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We study the decays of ultraheavy [nix > 10^"^ GeV) and quasistable (lifetime tx much larger 
than the age of the Universe to) particles as the source of Ultra High Energy Cosmic Rays (UHE 
CR). These particles are assumed to constitute a tiny fraction ^x of CDM in the Universe, with 
being the same in the halo of our Galaxy and in the intergalactic space. The elementary-particle and 
cosmological scenarios for these particles are briefly outlined. The UHE CR fluxes produced at the 
decays of X- particles are calculated. The dominant contribution is given by fluxes of photons and 
nucleons from the halo of our Galaxy and thus they do not exibit the GZK cutoff. The extragalactic 
components of UHE CR are suppressed by the smaller extragalactic density of X-particles and hence 
the cascade limit is relaxed. We discuss the spectrum of produced Extensive Air Showers (EAS) 
and a signal from Virgo cluster as signatures of this model. 
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The observations of Ultra-High Energy Cosmic Rays 
(UHE CR) reveal the presence of a new, isotropic com- 
ponent at energies E > 1- 10^*^ GeV (for a review see Ref. 

). This component is thought to have an extragalactic 
origin since the galactic magnetic field cannot isotropize 
the particles of such energies produced by astrophysical 
sources in the Galaxy. On the other hand, the observa- 
tion of particles of the highest energies, especially of the 
two events with energies 2 — 3 • 10" GeV g, contradicts 
the GZK cutoff §1 at £; ~ 3 ■ 10^° GeV, which is the sig- 
nature of extragalactic UHE CR. All known extragalactic 
sources of UHE CR, such as AGN [Q , topological defects 
1^ or the Local Supercluster result in a well pro- 
nounced GZK cutoff, although in some cases the cutoff 
energy is shifted closer to 1 • 10^^ GeV ||]. UHE neu- 
trinos could give a spectrum without cutoff, but the 
neutrino fluxes and the neutrino-nucleon cross-section are 
not large enough to render the neutrino a realistic can- 
didate for the UHE CR events. 

In this Letter, we propose a scenario in which the UHE 
CR spectrum has no GZK cutoff and is nearly isotropic. 
Our main assumption is that Cold Dark Matter (CDM) 
has a small admixture of long-lived supermassive X- 
particles. Since, apart from very small scales, fluctua- 
tions grow identically in all components of CDM, the 
fraction of X-particles, ^x, is expected to be the same in 
all structures. In particular, is the same in the halo of 
our Galaxy and in the extragalactic space. Thus the halo 
density of X-particles is enhanced in comparison with the 
extragalactic density. The decays of these particles pro- 
duce UHE CR, whose flux is dominated by the halo com- 
ponent, and therefore has no GZK cutoff. Moreover, the 
potentially dangerous cascade radiation |^ is suppressed. 
Long-lived massive relic particles were already discussed 
in the literature as a source of high energy neutrino ra- 
diation However, in our case the particles must be 
much heavier {mx ~ 10^^ — 10^^ GeV). 

The plan of our paper is as follows. First, we take an 



phenomenological approach and treat the density nx of 
AT-particles and their lifetime tx as free parameters fixed 
only by the requirement that the observed UHE CR flux 
is reproduced. We calculate the fluxes of nucleons, pho- 
tons and neutrinos, considering the production of cascade 
radiation, positrons, antiprotons and radio fluxes as con- 
straints. We then discuss how the required properties of 
AT-particles can be realized. 

The decays of X-particles result in the production of 
nucleons with a spectrum W]s[{mx , x), where mx is the 
mass of the X-particle and x = E/mx- The flux of nu- 
cleons {p,p,n,n) from the halo and extragalactic space 
can be calculated as 

pN{E)^^'^R^—WNimx,x), (1) 
47r Tx rnx 

where index i runs through h (halo) and ex (extragalac- 
tic), Ri is the size of the halo Rh, or the attenuation 
length of UHE protons due to their collisions with mi- 
crowave photons, Xp{E), for the halo case and extragalac- 
tic case, respectively. We shall assume mxn\ = ■C^Pcdm 
and rnxn^x = ^x^cuMPcr, where ^x describes the frac- 
tion of X-particles in CDM, r^cDM is the CDM density 
in units of the critical density pcr- We shall use the fol- 
lowing values for these parameters: a large DM halo with 
Rh ~ 100 kpc (a smaller halo with R^ = 50 kpc is possi- 
ble, too), ^IcDMh'^ = 0.2, the mass of X-particle in the 
range 10^^ QgV < mx < lO^^ GeV, the fraction of X- 
particles ^ 1 and tx 3> toj where to is the age of 
the Universe. The two last parameters are convolved in 
the flux calculations in a single parameter rx — S,xto/Tx- 
Following [0 , we shall use the QCD fragmentation func- 
tion in MLLA approximation (see [pd|) 

W7v(TOx,a;) = — exp ( ^ — 1, (2) 
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where 
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X = E/mx, = (A/mx)^/^ and A = 0.234 GeV. The 
normahzation constant Kn is found from energy conser- 
vation as 



Kn I dx exp 
'o 



In^ x/xy, 
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where /jv is the fraction of energy transferred to nucle- 
ons. Using Z°-decay as a guide, we assume Jn ~ O.OS/tt, 
where is the corresponding fraction for pions (LEP 
gives 0.027 for pp only). For the attenuation length of 
UHE protons due to their interactions with microwave 
photons, we use the values given in the book 

The high energy photon flux is produced mainly due 
to decays of neutral pions and can be calculated for the 
halo case as 



47r Tx 



(3) 



where Nj(E) is the number of photons with energy E 
produced per decay of one X-particle. The latter is given 

by 



N^{E) 
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The normalization constant K^^o is again found from the 
condition that neutral pions take away the fraction /tt/S 
of the total energy mx- 

For the calculation of the extragalactic gamma-ray 
flux, it is enough to replace the size of the halo, Rh, 
by the absorption length of a photon, A-y (E) . The main 
photon absorption process is e+e^ pair production on 
background radiation and, at i? > 1 • 10^° GeV, on the 
radio background. The neutrino flux calculation is simi- 
lar. 

Before discussing the obtained results, we consider var- 
ious astrophysical constraints. 

The most stringent constraint comes from electromag- 
netic cascade radiation, which is initiated by high-energy 
photons and electrons from pion decays and is developing 
due to interaction with low energy background photons. 
The relevant calculations were performed in Ref. In 
our case this constraint is weaker, because the low-energy 
extragalactic nucleon flux is ~ 4 times smaller than that 
one from the Galactic halo (see Fig. 1). Thus the cascade 
radiation is suppressed by the same factor. 

The relevant parameter which characterizes the flux of 
cascade radiation is the total energy density of cascade 
radiation Wcas- The observation of the low-energy diffuse 
gamma-ray flux results in the limit Wcas < 1 • 10~^ — 1 ■ 
10~^ eV/cm^ B. In our case, the cascade energy density 



calculated by integration over cosmological epochs (with 
the dominant contribution given by the present epoch 
z — 0) yields 
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6.3 • lOVx/x eV/cm^ (5) 



To flt the UHE CR observational data by nucleons from 
halo, we need rx = 5 ■ 10"^-^. Thus the cascade energy 
density is Wcas = 3.2 • 10~^/7r eV/cm'^, well below the 
observational bound. 

The other constraints come from the observed fluxes 
of positrons and antiprotons in our Galaxy and from the 
isotropic component of the radio flux. We performed 
detailed calculations which will be published elsewhere. 
In all cases the abovementioned constraints are satisfled 
and they are weaker than that due to cascade gamma- 
radiation. 

Now we address the elementary-particle and cosmolog- 
ical aspects of a supermassive, long-living particle. Can 
the relic density of superheavy X-particles be as high as 
required in our calculations? And can this particle have 
a lifetime comparable or larger than the age of the Uni- 
verse? 

Let us assume that X is a neutral fermion which be- 
longs to a representation of the SU{2) x U{1) group. We 
assume also that the stability of AT-particles is protected 
by a discrete symmetry which is respected by all interac- 
tions except quantum gravity through wormhole effects. 
In other words, our particle is very similar to a very heavy 
neutralino with a conserved quantum number, i?', being 
the direct analogue of i?-parity (see Q and the refer- 
ences therein). Thus, one can assume that the decay 
of X-particles occurs due to dimension 5 operators, in- 
versely proportional to the Planck mass mpi and addi- 
tionally suppressed by a factor exp(— 5'), where S is the 
action of a wormhole which absorbs one i?'-charge. As 
an example one can consider a term 
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where 5* describes X-particle, and is a SU{2) scalar 
with vacuum expectation value vew = 250 GeV. After 
spontaneous symmetry breaking the term (^ results in 
the mixing of A"-particle and neutrino, and the lifetime 
due to X—>-iy + q + q, e.g., is given by 
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192(27r) 
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where Gi? is the Fermi constant. The lifetime tx > to 
for X-particlc with mx > 10^^ GeV needs S > 44. This 
value is within the range of the allowed values as dis- 
cussed in Ref. . 

Let us now turn to the cosmological production of X- 
particles with mx > 10^^ GeV. Several mechanisms can 
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be considered, including thermal production at the re- 
heating stage, production through the decay of inflaton 
field at the end of the "pre-heating" period following in- 
flation, and through the decay of hybrid topological de- 
fects, such as monopoles connected by strings or walls 
bounded by strings. 

For the thermal production, temperatures comparable 
to mx are needed. In the case of a heavy decaying grav- 
itino, the reheating temperature Tfj (which is the highest 
temperature relevant for our problem) is severely limited 
to value below 10^ — 10^° GeV, depending on the gravitino 
mass (see Ref. and references therein) . On the other 
hand, in models with dynamically broken supersymme- 
try, the lightest supersymmetric particle is the gravitino. 
Gravitinos with mass 7713/2 < 1 keV interact relatively 
strongly with the thermal bath, thus decoupling rela- 
tively late, and can be the CDM particle In this sce- 
nario all phenomenological constraints on Tj^ (including 
the decay of the second lightest supersymmetric particle) 
disappear and one can assume Tr ^ 10^^ — 10^^ GeV. In 
this range of temperatures, X-particles are not in thermal 
equilibrium. If Tr < mx , the density nx of X-particles 
produced during the reheating phase at time Ir due to 
a + d^X + X\s easily estimated as 

nxitR) ^ Naulexp^^-"^:^^^ (TxtR, (8) 

where Na is the number of flavors which participate in the 
production of X-particlcs, is the density of a-particles 
and ax is the production cross-section. The density of 
X-particles at the present epoch can be found by the 
standard procedure of calculating the ratio nx/ s, where 
s is the entropy density. Then for mx = 1 • 10^^ GeV and 
£,x in the wide range of values 10~^ — 10^^, the required 
reheating temperature is Tr 3 ■ 10^^ GeV. 

In the second scenario mentioned above, non- 
equilibrium inflaton decay, X-particles are usually over- 
produced and a second period of inflation is needed to 
suppress their density. 

Finally, X-particles could be produced by the decay of 
hybrid topological defects, i.e. monopoles connected by 
strings or walls bounded by strings. For example, strings 
of energy scale rjs ^ tiix could be formed at a phase tran- 
sition at or near the end of inflation. At a second phase 
transition with symmetry-breaking scale rjw < "mx each 
string gets attached to a domain wall. The wall ten- 
sion pulls the strings together and eventually leads to 
a breakup of the network. The resulting pieces of wall 
bounded by string lose their energy by gravitational ra- 
diation and by particle production. X-particles are pro- 
duced whenever strings cross one another and also in the 
decay of the pieces which fragmented down to the size 
comparable to the string thickness. The X-particle den- 
sity produced in this way depends on the details of the 
fragmentation process, but rough estimates suggest that 



the required values of nx / s can be obtained for a wide 
range of string and wall parameters. 

Let us now discuss the obtained results. The fluxes 
shown in Fig. 1 are obtained for Rh — 100 kpc, mx = 
1-1013 GeV and rx = ixh/rx = S-IO"". This ratio rx 
allows very small and tx > to. The fluxes near the 
maximum energy -Emax = 5 • 10^^ GeV were only roughly 
estimated (dotted lines on the graph). 

It is easy to verify that the extragalactic nucleon flux 
at i? < 3 • 10^ GeV is suppressed by a factor ~ 4 and 
by a much larger factor at higher energies due to nucleon 
energy losses. The flux of extragalactic photons is sup- 
pressed even stronger, because the attenuation length for 
photons (due to absorption on radio-radiation) is much 
smaller than for nucleons (see Ref. fl^]). This flux is not 
shown in the graph. The flux of high energy gamma- 
radiation from the halo is by a factor 7 higher than that 
of nucleons and the neutrino flux, given in the Fig.l as the 
sum of the dominant halo component and subdominant 
extragalactic one, is twice higher than the gamma-ray 
flux. 

The spectrum of the observed EAS is formed due to 
fluxes of gamma-rays and nucleons. The gamma-ray 
contribution to this spectrum is rather complicated. In 
contrast to low energies, the photon-induced showers at 
E > 10^ GeV have the low-energy muon component as 
abundant as that for nucleon-induced showers . How- 
ever, the shower production by the photons is suppressed 
by the LPM effect and by absorption in geomagnetic 
field (for recent calculations and discussion see [^,0 and 
references therein). These effects are energy dependent. 
The LPM effect starts at 10^ - lO^" GeV and it al- 
most fully suppresses the production of " normal" EAS at 
> 1 ■ 10^^ GeV, when maximum EAS reaches the see 
level practically for all zenith angles The calculation 
of the spectrum of EAS is outside the scope of this paper 
and the normalization of the halo nucleon spectrum by 
observational data at i? ~ 2 • 10^^ GeV in Fig. 1 has an 
illustrative character. The general tendency of greater 
suppression of photon-induced showers with increase of 
energy might improve the agreement between calculated 
and observed spectra. 

We wish to note that the excess of the gamma-ray flux 
over the nucleon flux from the halo is an unavoidable 
feature of our model. It follows from the more effective 
production of pions than nucleons in the QCD cascades 
from the decay of AT-particle. 

Although AT-particles with necessary properties can be 
produced by a variety of mechanisms, it should be noted 
that their lifetime and spatial density had to be fine- 
tuned to the desired values with the help of exponential 
factors. 

The signature of our model might be the signal from 
the Virgo cluster. The virial mass of the Virgo clus- 
ter is Afviigo ~ 1 ■ IO^^Mq and the distance to it 
ii = 20 Mpc. If UHE protons (and antiprotons) prop- 
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agate rectilinearly from this source (which could be the 



case for En 
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10 GeV), their flux is given by 



^Virgo ^ 



'^WN{mx,x) 



toR^r 



(9) 



The ratio of this flux to the diffuse flux from the half 
hemisphere is 6.4 • 10~^. This signature becomes less 
pronounced at smaller energies, when protons can be 
strongly deflected by intergalactic magnetic fields. 

When our work was in progress, we learned that a 
similar idea was put forward by V. A. Kuzmin and 
V. A. Rubakov |2^. The main difference is that we 
take into account the radiation from the galactic halo, 
which is the main issue of our work, while the authors 
above limited their consideration to the extragalactic 
component. We are grateful to V.A. Kuzmin and V.A. 
Rubakov for interesting discussions. M.K. was supported 
by a Feodor-Lynen scholarship of the Alexander von 
Humboldt-Stiftung. 
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FIG. 1: Predicted fluxes from decaying X-particles: nucleons {p,p, n, n) from 
the halo (curve extragalactic protons (curve /p''*'^), photons from the halo 

(curve /^*'°), and neutrinos from the halo and the extragalactic space (curve 
jtot-j rjjjg data points are based on the compilation made in Ref. [21]. 



